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Abstract After a brief reminder of the mechanism of gravitational lensing of extended 
sources, the particular case of the host galaxy of QSO RXJ091 1, a high redshift (z ~ 2.8) 
quadruply imaged quasar, is explored. The non linearity of the problem, together with the 
proximity of the source to a cusp of the lens inner caustic, have important consequences 
on the dependence of the image appearance on the size and shape of the source. Their 
expected main features and their interpretation in terms of source extension and shape are 
investigated in a spirit of simplicity and in preparation for the analysis of high sensitivity 
and spatial resolution images that will soon be within reach with the completion of the 
Atacama Large Millimeter/submillimeter Array (ALMA). In particular, the information 
on source size carried by relative image brightness is discussed. Extension of the results 
to other types of quadruply imaged quasars is briefly considered. 
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1 INTRODUCTION 

Strong gravitational lensing has become a textbook topic. In particular, several authors, such as 
Blandford et al. d!989t or Saha and Williams (120031 . have summarized the main properties in simple 
terms, underlining the most general qualitative features. While the case of complex lens configurations 
has been extensively studied (see for example Julio et al. ( 120071 1), in particular with the aim of evaluating 
the mass distribution of baryonic and dark matter in cluster lenses, studies of strongly lensed extended 
sources are less common. The clearest cases of strong lensing, which are naturally the most studied, 
are often associated with sources located near the inner caustic of the lens, making the problem highly 
non linear: when crossing the caustic outward, magnifications become infinite and one switches from a 
four-image to a two-image configuration. Several authors, such as Dominik (1988), Bartelmann ( 120031 1, 
Suyu and Blandford (120051 1 or Suyu et al. (2009) have analysed the consequences in the case of ex- 
tended sources and described their effects in some detail. With the completion of high sensitivity and 
angular resolution instruments, such as the Atacama Large Millimeter/submillimeter Array (ALMA) or 
the Square Kilometer Array (SKA), the study of host galaxies of high redshift gravitationally lensed 
quasars will enter a new phase with the ability of resolving spatially the sources in both their molecular 
gas and dust contents. 

The present work concentrates on the case of a quadruply imaged quasar, RXJ0911, which pro- 
vides a good illustration of the main features. Observations in the visible and near-infrared by Burud et 
al. JT998), using the 2.56 m Nordic Optical Telescope and the ESO 3.5 m New Technology Telescope, 
resolved clearly the object in four QSO images having a redshift z ~ 2.8 and an elongated lens galaxy 
and measured the relative image fluxes, their magnitudes and positions. Three of the four images, named 
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CASTLES 



Fig. 1 Image configuration (Kochanekl999j of RXJ091 1 showing the lensing galaxies G and 
G'. 



Table 1 Position data from CASTLES (upper line, Kochanek et al. 11 9991 and best fit results 
from Burud et al. ( 1998 ). lower line. Magnifications (M) are from Burud et al. (119981 model. 



Observations 


Al 


A2 


A3 


B 


G 


G' 


Position 


-RA(arcsec) 






-0.260 
-0.259 


0.018 
0.013 


2.972 
2.935 


0.698 


1.452 


Dec(arcsec) 






0.406 
0.402 


0.960 
0.946 


0.792 
0.785 


0.512 


1.177 


M 




-4.45 


8.59 


-3.70 


1.79 







Al to A3, are at similar distances from the lens (~ 0.9") while the fourth image, B, is 2.5 times as far 
(~ 2.2"). Modelling such lensing properties (Burud et al. 119981 requires the presence of an external 
shear in addition to the lens ellipticity. Two years later, the lens was identified by Kneib et al. ( 120001 ) as 
being a member of a galaxy cluster having a mean redshift z = 0.769 ± 0.002 and a velocity dispersion 
of 836 ± 190 km s _1 using the Low Resolution Imaging Spectrograph of the Keck II Telescope. The 
cluster has an estimated mass of ~ (6.2 ± 2.7)10 14 solar masses. 

Recently, RXJ0911 has been observed using the Plateau de Bure Interferometer (Tuan Anh et 
al. 120121 1 tuned in the CO(7-6) region. At such wavelengths, the source is no longer the quasar, which 
can be considered as a point source, but its host galaxy. More precisely, the CO line is associated with 
the presence of molecular gas and the continuum with that of dust: measuring the source extension in 
each case would provide important information on the physics of early galaxies. The Plateau de Bure 
data resolve the source spatially on the CO line, suggesting a root mean square radius of several hundred 
parsecs (1" in the sky corresponds to 8 kpc in the source plane). The evidence for resolving spatially 
the source in the continuum is much weaker, at the level of only two standard deviations. Such results 
exploit the performance of the radio telescope array to the limits of its ability but give confidence that 
the better acceptance, resolution and sensitivity of ALMA, soon to be completed, will give more defini- 
tive results. The aim of the present work is to describe in simple terms what can be expected in order to 
guide the analysis of future higher resolution data. 

A HST image of the quasar is shown in Figure Q] the associated data are listed in Table 1. 

2 STRONG LENSING: A REMINDER 
2.1 General formalism 

Fermat principle states that images form where the gradient of the time delay, r, cancels. In the approx- 
imation of small deflections, which always applies in practice, the time delay can be written as the sum 
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of a geometrical delay and of the gravitational delay proper: r = ro[l/2(i — s) 2 — ip]. Here, i ands are 
the image and source vectors in sky coordinates (in a plane normal to the line of sight), To is a constant 
time scale and ip is an effective potential that describes the deflection induced by the lens as a function 
of the sky coordinates of the image. The effective potential ip is proportional to the integral of the gravity 
potential along the line of sight between source and observer. A convenient form, used by many authors, 
includes an elliptical lens and an external shear, the axes of the ellipse being taken as coordinate axes 
without loss of generality: 

ip = r r(l + ecos2ip) 1/2 + l/27 r 2 cos2(y> - ip ) , (1) 

where (r, ip) are the polar coordinates of the image. The potential is usually given in Cartesian coordi- 
nates, as in Burud et al. (1998 ); we rewrote it here in polar coordinates, which we found convenient in 
the case where the potential takes such a simple form, but this is a technical detail of no relevance to 
the results of the study. The lens term, of strength tq and aspect ratio [(1 + e)/(l — e)] 1 ^ 2 , decreases 
as 1/r outside the core region. The shear term has a strength 70 and makes an angle tpo with the major 
axis of the lens ellipse. Writing that the gradient of Relation 1 cancels, and calling (r s , ip s ) the polar 
coordinates of the source, one obtains the lens equation: 

r s e llps = re ltp (l - r^dip/dr - ir~ 2 dip/dip) , (2) 



As is well known, there may typically be two or four images depending on the position of the source 
with respect to the inner caustic of the lens. If the potential is isotropic, the lens equation reduces 
to r s e lips = e lip (r — dip/dr), which has two obvious solutions, one at ip + = ip a and the other at 
(p_ = (p s + 7r with r + = dip/dr + r s and r_ = dip/dr — r s respectively. For r s = 0, the alignment is 
perfect and one obtains an Einstein ring having r = dtp /dr. 



2.2 Extended sources 



To the extent that the source is small and not too close to the lens inner caustic, the image of an extended 
source is simply obtained by differentiating the lens equation, r s e lips = Dre ltp , where D = D r + iDi 
and D r = 1 — r~ 1 dip / 'dr, D L = —r~ 2 dip / dp . We obtain this way the relation between a point 
(r s + dr s , ip s + dtps) on the source and its image (r + dr, cp + dip): 

(dr s + ir s d(p s )e llps = D(dr + irdtp)e lip + (^-dr + ^-dip)re iip (3) 

dr dtp 

In practice, one does not directly observe the source, but only its images. It is therefore convenient to 
replace, in the left-hand side of the above relation, the source dependent term r s e lvs by its expression 
in terms of the image coordinates: 

D(^+idcp s ) = (D + r^)- + (iD + ^-)d V (4) 
r s dr r dip 

Relation 3 gives the coordinates of an image point as a function of those of the corresponding source 
point. Indeed, dr s and r s dip s are Cartesian coordinates having their origin at the centre of the source and 
the axis of abscissas radial outwards; similarly, dr and rd<p are Cartesian coordinates having their origin 
at the centre of the image and the axis of abscissas radial outwards (Figure |2}. For Relation 3 to apply, 
dr s and r s d<p s must be small enough for the corresponding source points to stay away from the lens 
inner caustic. In practical cases, with the centre of the source located near the caustic, this will often 
not be the case if the source extension is such that it overlaps the caustic. The linear approximation 
of Relation 3 needs therefore to be used with care. Yet, it usefully serves several purposes, such as 
providing explicit expressions for the magnifications which may be calculated using arbitrarily small 
values of dr s and r s d(p s - or for giving a qualitative illustration of the main features in simple terms 
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Fig. 2 Left panel: Schematic geometry and definition of coordinates. The origin is the centre 
of lens G, the polar axis is positive westward and angles are measured counter clockwise. 
Right panel: Isotropic potential imaging. 



as is done below. Relation 3 being linear, it is straightforward to express (dr/r, dip) as a function of 
(dr s /r s ,d<p s ) in terms of a lensing matrix A defined as: 

dr/r = \\\dr s /r s + \i 2 dip s (5a) 
dp = X 2 idr s /r s + X 22 dp s (5b) 

with: 

A u = [D r {D r + dDi/dip) + A(A - dD r /d<p)]/n 
A12 = -(D r dD r /d<p + DidDi/d<p)/n 
A 2 i = r(DidD r /dr - D r dDi/dr})/n 

A 22 = [D r (D r + rdD r /dr) + D t (D z + rdDi/dr)]/ f i W 
(j, = (D r + rdD r /dr){D r + dDi/dip)+ 
(A + rdDi/dr){Di - dD r /d<p) 

For an isotropic source brightness, such that < dr s /r s >= 0, < dip s >= 0, < dp s dr s /r s >= 
and < dr 2 s >= r 2 < dip 2 >= a 2 , it is straightforward to write the expressions of the semi-major and 
semi-minor axes a and b of the image, its position angle 8 and magnification M: 

{a±b) = (ar/r s )[|A| 2 ± 2det{\)] 1 / 2 

tan20 = (A11A21 + A 12 A 22 )/(A? 1 + \\ 2 - \ 2 21 - A| 2 ) (7) 
M = (r/r s ) 2 det{\) , 

where det(X) = AnA 22 — A i2 A 2 i and |A| 2 = Afj + A 2 2 + A 21 + X 22 . 

The case of an isotropic potential (Figure [2] right) illustrates the above relations, with tan26 = 0, 
a = a(l — d 2 t/j/dr 2 )~ 1 and b = ar/r s . The image is stretched normally to the lens-image direction 
by a factor r/r s . The magnification in the lens-image direction is unity when the second derivative of 
the potential cancels. When small perturbations are added to the isotropic potential in the form of an 
external shear and/or a non-zero eccentricity of the lens, a, b and 8 deviate by correspondingly small 
amounts from the above values. When r s decreases, the two images get more and more elongated in the 
tangential direction forming characteristic arcs and finally merge while drifting toward the Einstein ring 
of radius dip /dr. 

The case of a potential ip = rro(l + £cos2</3) 1 / 2 that describes an elliptical lens without external 
shear, e accounting for the ellipticity of the lens, is particularly simple: one then has < dr 2 >= a 2 , 
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Fig. 3 Dependence of the image positions and morphology on the source position. The lens 
is circular (rrj = 1") and centred at the origin. The potential includes an external shear, 
70 = 0.2, if = 0. From left to right: tp s = 0°, r s = 0.09", 0.18", 0.27", 0.36", 0.45" 
(upper panels); ip s = 45°, r s = 0.06", 0.12", 0.18", 0.24", 0.30" (second row); tp s = 90°, 
r, = 0.10", 0.25", 0.40", 0.55", 0.70" (lower panels). The source is a disk of radius 0.02". 

< r 2 dip 2 >= er 2 /yu 2 and < rdrdtp >= 0. The images are ellipses stretched tangentially, their minor 
axes are equal to the source diameter and, writing fj, = 1 — ror _1 (l — e 2 )(l + £cos2</?)~ 3 / 2 , their major 
axes are a factor 1/fj, larger. At constant <p, 1/fj, is of the form r/(r — kro) with k = (1 — e 2 )(l + 
ecos2<^) -3 / 2 : it reaches unity at large values of r and increases when r decreases toward the critical 
curve, of equation r = kro, where l//x diverges. 

Similar results are obtained for a potential of the form -0 = t-qt + 1 / 2jor 2 cos2ip that describes 
a circular lens with external shear along the axis of abscissas. Figure [3] shows how the positions of 
the images vary as a function of the source position: one switches from a four-image configuration, 
when the source is inside the caustic, to a two-image configuration, when the source is outside the 
caustic. Image shapes, illustrated in Figure |4]for different source positions, are again ellipses stretched 
tangentially but, at variance with the preceding case, they are slightly tilted as soon as tp s deviates from 
0° or 90°. Qualitatively, for small values of the external shear, the general features of the isotropic case 
are maintained. 

3 QSORXJ0911 

3.1 Quasar point images 

We study the image configuration of RXJ091 1 using one of two potentials: 

V>i = r r[l + ecos2(<^ - <^i)] 1/2 + l/27 r 2 cos2(</? - ip Q ) (8) 

use d by m any authors (Blandf ord & Kochanek [19871 Sa ha & Willia ms [20031 Peeples [20031 Burud et 
al. ([19981 . Kassiola & Kovner[[993l Witt & Mao [T997l Kneib et al. [2000b and including an elliptical 
main lens G and an external shear, or 

i'2 = r r + r' [r* 2 + r 2 - 2rr*cos(tp - ip*)] l/2 + l/2j r 2 cos2(tp - tp ) (9) 
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Fig. 4 Image appearances for the same potential and source size as in Figure [3] Source posi- 
tions are: <p a = 20°, r s = 0.1", 0.15", 0.2", 0.25", 0.3" from left to right. Images are in [dr 
(up), rdip (left)] coordinates. Images 1 to 4, numbered from South clockwise, are displayed 
in the upper to lower rows respectively. 



describing a circular main lens G, a small perturbation from its satellite galaxy G' and an external shear 
[KochanekHHH SchechterUxBl- 

The parameters of the model are adjusted by minimizing the x 2 describing the match between the 
observed quadruple point images and the prediction of the point source model. The reason for using two 
different potentials is to illustrate the relative independence of most of the main results on the detailed 
form of the potential, not to argue about which form is better. Care has been taken, when exploring the 
X 2 dependence on the coordinates of the point source, to search for minima confined within the lens 
inner caustic. This precaution is important when the minimization algorithm approximates locally the 
surface of constant x 2 by a quadratic form in order to evaluate the gradient along which to move and to 
calculate uncertainties. The results are summarized in Tables 2 and 3. Images obtained using potential 
ipi for a disk source of uniform brightness over a radius R s = 0.025" are shown in Figure|5] together 
with the critical curve. The images obtained using potential ip2 are practically indistinguishable from 
those obtained using potential ipi . 

3.2 Extended source 

As can be seen from Table 1, the distance between the centres of image A2 and of its neighbours is 
~ 0.6" and the sums of the semi-major axes of A2 and Al or A3 is ~ 9 source radii: as soon as 
the source radius exceeds ~ 0.06" to 0.07", the three A images merge and the linear approximation 
developed in the preceding section is no longer appropriate. In the present section, images are therefore 
constructed point-by-point, each image point being given a weight equal to the absolute value of its 
magnification. In practice, we calculated, for each point of a fine-mesh map of the source plane, the 
positions and magnifications of its two or four image points in the image plane and we constructed 
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Table 2 Best fit results (r in arcseconds and ip in degrees). 





Al 


A2 


A3 


B 


fobs 


0.866 


0.964 


0.814 


2.291 


r m odel{lpl) 


0.853 


0.968 


0.819 


2.294 


r m odel[lp2) 


0.845 


0.967 


0.828 


2.282 


<Pobs 


216.3 


186.3 


146.6 


7.0 


•Pmodel (i>l) 


219.3 


188.2 


150.0 


5.8 


•Pmodel (i>2) 


219.0 


188.4 


151.2 


6.2 


M(^l) 


-4.14 


8.33 


-3.37 


1.92 


M(V> 2 ) 


-4.45 


8.32 


-3.47 


1.76 



Table 3 Best fit values of the parameters (angles in degrees and angular distances in arcsec- 
onds). 



Potential ipi 


ro 




<pi 


7o 






(fis 


1.1085 


-0.0237 


65.0 


0.309 


7.32 


0.4456 


3.84 


Potential ip 2 (r* = 1.005, tp* = 41.4) 


ro 






70 


tpo 






1.0931 


0.01935 




0.317 


6.7 


0.4500 


6.7 
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Fig. 5 Image positions obtained using potential ifti are shown in blue for a source radius of 
0.025". The source is shown in red. The lens G is at the origin of coordinates (black asterisk). 
Red crosses indicate the HST observed image positions. The results obtained using potential 
ip2 are identical within the precision of the figure. Sky coordinates are in arcsec. 



images of a set of source points as the set of images of these source points, each time using the proper 
weight (including magnification and, in case of a non-uniform distribution, the source brightness). Both 
potentials give very similar results (Figure |6}, with the source covering the western cusp of the caustic: 
the inner part of it gives four images, the outer part, only two. Indeed, as commented earlier by Saha 
& Williams ( 120031 1. the precise position of the point images, including the strong asymmetry between 
the A and B images, severely constrains the position of the source with respect to the caustic cusp, with 
little dependence on the precise form of the effective potential. Uniform (radius R s ) and G aussian (<r s ) 
brightness distributions are also found to give similar results in terms of p s = \J < r 2 s > = R s /y/2 = 
a s y/2. Figure|7]displays the images obtained using potential ip\ for different source extensions. 

The relative occurrence of two-image configurations increases with the size of the source, causing 
the A images to become globally fainter with respect to image B, the latter being always present in the 
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Fig. 6 Position of the source with respect to the caustic for potentials ipi (left panel) and ip2 
(right panel). The concentric circles are for R s — 0.05", 0.1", 0.2" and 0.3". The white areas 
show the inner region of the caustic where the source needs to be located to form four images. 
When the source is outside this region only two images are formed. Sky coordinates are in 
arcsec. 
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Fig. 7 Images of RXJ0911 (assuming a disk source of radius R s = 0.05", 0.1", 0.15" and 
0.2" from left to right). The white area shows the caustic and the source is shown as a red 
disk. Sky coordinates are in arcsec. 



two-image configuration. This suggests using the apparent integrated flux ratio, B/A, as a measure of 
the source size. This ratio is essentially unaffected by beam convolution and its measurement does not 
require a very good angular resolution (yet sufficient to separate A from B). However, as image B is not 
much magnified, it requires a high sensitivity to be obtained with good precision. 

To understand what happens, it is important to keep in mind that magnifications depend strongly on 
the position of the source point with respect to the caustic cusp, or equivalently on the position of the 
image point with respect to the critical curve. This is illustrated in Figure|8] which displays the variations 
over the source plane of the A and B magnifications separately. 

The A images (considered here globally, namely the flux being integrated over the area covered 
by the three A images) are faint in most of the two-image region, large magnifications being only 
reached in the vicinity of the caustic cusp. On the contrary, the magnification of image B, which is 
far from the critical curve, is always weak, whether the source point is inside or outside the caustic 
cusp: it varies very slowly across the region explored in Figure [8] When the source size increases, 
starting from a B/A integrated flux ratio of ^1 1% for a point source, as soon as the image size overlaps 
sufficiently the caustic, the B integrated flux grows essentially in proportion with the source area, the 
average magnification remaining between 1.5 and 2.2. But the relative growth of the A integrated flux 
is much slower: outside the caustic cusp, there is only one A image left and the magnification decreases 
rapidly when moving away from the cusp. Hence an increase of the B/A integrated flux ratio as a 
function of source size as displayed in the right panel of Figure [8] 



Gravitationally lensed extended sources: the case of QSO RXJ091 1 



9 




Fig. 8 Distribution of the magnifications in the source plane for the three A images considered 
globally (left panel) and for image B (middle panel). Concentric circles show sources of 
respective radii R s = 0.01", 0.06", 0.11", 0.16", 0.21", 0.26" and 0.31". In the left panel, 
from left to right, contours are for M = 6, 8, 10 and 12. In the middle panel, again from left to 
right, contours are for M = 2.0, 1.9, 1.8, 1.7 and 1.6. The right panel displays the dependence 
of the B/A integrated flux ratio (ordinate) on p s = R s /^/2 (arcseconds) for a disk source of 
uniform brightness. One arcsecond corresponds to 8 kpc in the source plane. 

3.3 Different configurations 

The above comments apply, mutatis mutandis, to other cases of quadruply imaged quasars. Using the 
classification proposed by Saha & Williams (120031 1 we briefly consider a few typical configurations: 
a "core quad", such as Q2237+030, H1413+117 (the Cloverleaf) or HST 1411+521, with the source 
close to the origin; an "inclined quad", such as MG0414+0534 or B1608+656, with the source close 
to the caustic but far from its cusps; a "long axis quad", such as B 1422+231, with the source near a 
major axis cusp of the caustic. The case of RXJ091 1, studied here in detail, is a "short axis quad" with 
the source near the minor axis cusp. For each of these typical configurations, we use potential ipi with 
the parameters taking the same values as for RXJ091 1 (Table 2) and simply change the position of the 
source. 

The results are displayed in Figure [9] In the case of a long axis quad, as in the present case of a 
short axis quad, the B/A integrated flux ratio provides a good measure of the source size, however with 
a much lower magnification of image B, resulting from its being far from the critical curve: in such a 
long axis case, both A and B magnifications depend strongly on the distance of the source to the cusp. 
In the short axis case of Figure [9] this distance is 0.11" while in the long axis case, two values, 0.07" 
and 0.35" are displayed. In the case of a core quad, the source does not overlap the caustic, even at the 
maximal value of p s and the images do not merge: their relative brightness is almost independent from 
the source size but their radial extensions provide a good measure of it. Diametrically opposite images 
have nearly equal magnifications that are averaged in Figure [9] 

In the inclined quad case, the two north-eastern images merge as soon as the source overlaps the 
caustic and the ratio of the global brightness of the south-western images to that of the north-eastern 
images in displayed in Figure [9] From this analysis we conclude that integrated flux ratios provide a 
useful measure of the source size in configurations where it is located in the vicinity of the inner caustic. 
However, when it is far from it, in configurations that are closer to an Einstein ring as is the case for the 
Cloverleaf, image extensions must be preferred to integrated flux ratios to evaluate the source size. 

3.4 Velocity gradients 

Finally, we briefly comment on the effect of a velocity gradient on the appearance of the A images. Such 
gradient can be expected, for example, from a movement of rotation of the source or from molecular 
outflow. In the absence of gravitational lensing, the maps of the parts of the source associated with 
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Fig. 9 Typical quadruply imaged quasars according to Saha & Williams (2003;) classification. 
From left to right: short axis quad, long axis quad, core quad and inclined quad. The images 
(blue areas) are displayed in the upper panels for a disk source (red disk) having a radius 
R s = 0.05"; the caustic is shown as a grey area. The lower panels illustrate the dependence 
of image integrated flux on source size. The two left panels display the dependence of the B/A 
integrated flux ratio on p s with distances from the source centre to the cusp of 0.11" for the 
short axis and 0.35" (black) and 0.07" (red) for the long axis. The core quad panel displays 
the mean radial size of the images (NS in red and EW in black) as a function of p s . The right 
panel displays the brightness ratio of the south-western pair to the north-eastern pair. 



the red-shifted and blue-shifted parts of the velocity spectrum are obtained directly. However, in the 
presence of gravitational lensing, such is not the case when the source is extended. Reconstructing the 
maps of an extended source in general, and in particular of the parts of the source associated with 
the red-shifted and blue-shifted parts of the velocity spectrum, is more involved: one starts with some 
hypotheses about these maps, reconstructs the associated images and improves the hypotheses until 
one obtains a good match with the observed images. This iteration process is usually implemented by 
assuming some model for the source maps, the parameters of which are adjusted by minimizing the \ 2 
that measures the quality of the match between the observed and predicted images. At the end, if the 
models used for the blue-shifted and red-shifted parts of the source are adequate, the result is the same 
as in the absence of gravitational lensing. 

We illustrate this point in Figure [10] in the case of a disk source having R s = 0.15" and split 
in two halves by a diameter making different angles with the east-west direction. In the absence of 
gravitational lensing, each half disk of the source would have as image a similar half disk while in 
the presence of gravitational lensing, the images are distorted and do not provide a direct picture of 
the source. In particular, in the present case of three merging images, A2, which has a magnification 
about twice those of Al and A3, is inverted with respect to Al and A3, increasing the importance 
of distortion. This complication, inherent to the gravitational lensing of extended sources, is a minor 
difficulty in comparison with the important benefit offered by the large magnifications associated with 
gravitational lensing, without which the source could not even be resolved in most cases. 

4 SUMMARY 

The mapping of an extended source onto its gravitationally lensed images has been studied with the aim 
of guiding the analysis of future high resolution observations of the host galaxy of QSO RXJ091 1 and 
of similar high redshift quasars: the results are general enough to be of interest to the study of other 
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Fig. 10 A images of a split source having R s = 0.15". The split is along a diameter with 
position angles of 0°, 30°, 45°, 60°and 90°(from left to right). The upper panels are for the 
eastern half and the lower panels for the western half of the source. 

multiply imaged galaxies. Using simple effective potentials, the general features of the lensing mecha- 
nism have been recalled and explicit analytic relations have been given that are of help in understanding 
the main features. In particular Relations 5 and 6, describing the mapping in the linear approximation 
as a function of the lensing potential using a mapping matrix A, are useful in practical cases to calculate 
image magnifications and to provide a qualitative picture of the lensing mechanism. 

The study of the particular case of QSO RXJ0911 has illustrated how the problem becomes non 
linear when the source size becomes large enough for part of it to overlap the caustic. In such a case, 
the precise shape of the caustic in the source region becomes determinant in the construction of the 
images. The use of variables sensitive to the source size has been suggested, including the dependence 
of the relative image integrated fluxes. Typical quadruply imaged quasar configurations have been used 
as illustrations. The effect of strong lensing on a velocity gradient has been briefly considered. 
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